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Abstract 

Fermion self-energy associated with wave function renormalization is studied 
in a five-dimensional Yukawa theory on the orbifold S x jZi- One- loop divergence 
can be subtracted with only two renormalization constants in the bulk and on the 
branes. We show that the bulk and brane parts of the self-energy are uniquely 
determined with requiring physical conditions. With this procedure, momentum- 
scale dependence of the renormalized self-energy is given and the distribution of the 
bulk and brane parts can be found at low and high energies. Despite possible higher 
degrees of divergence in higher dimensions, the regularization scheme dependence 
does not arise. A viewpoint of the regularization scheme dependence at higher- 
loop level is also discussed. We find that the ratio of the bulk contribution to the 
brane contribution depends on the momentum scale in a very mild way, so that the 
relative coefficient of bulk and brane kinetic terms can be regarded as approximately 
constant for the leading quantum effect. The physical conditions given here are 
applicable to remove ambiguity in various orbifold models. 



1 Introduction 



The picture that physical quantities depend on the energy scale of interest is intuitive and 
has provided a clue to understand laws of Nature. A probe with a shorter wavelength in 
a system would resolve its substructure in a more microscopic way. If the substructures 
in the system are hierarchical, it can achieve predictability without knowledge of details 
at the other energy scales of the system. In a renormalizable quantum field theory in four 
dimensions, information for short and long distances is described well. Divergence that 
the theory can have for short distance is subtracted. After renormalization conditions are 
imposed, physical quantities are finite and they are energy-scale dependent. 

The energy-scale dependence of physical quantities is important not only for renor- 
malizable interactions but also for non-renormalizable interactions. Renormalizable in- 
teractions are treated with no new counterterms once all the renormalizable interactions 
are included. In usual four-dimensional models, non-renormalizable interactions are sup- 
posed to be suppressed by the ratio of the energy scale of interest to the ultraviolet 
momentum cutoff of the theory. In a four- dimensional theory where renormalizable and 
non-renormalizable interactions coexist, physical quantities can be dominated by contri- 
butions from lower-dimensional operators being renormalizable terms. In a theory with 
compactified extra dimensions, fields as four- dimensional modes can have dimension-four 
operators which are similar to renormalizable terms in four dimensions. From a standpoint 
that non-renormalizable interactions are irrelevant operators and that their contributions 
to physical quantities are negligible, it is interesting to search for rules or orders for 
possible effects in higher- dimensional field theory at each given loop level. 

In a theory with compactified extra dimensions, several characteristic properties need 
to be taken into account. They are higher- dimensional operators, regularization scheme 
dependence and brane terms. Even if the starting action integral includes all the operators 
up to a certain mass dimension, radiative corrections give rise to divergence for new local 
operators. Because higher-dimensional operators needed in the starting action integral 
affect the values of physical quantities, they must be determined in some ultraviolet 
completion or they must be found to be very small. As for scheme dependence, the 
point is that degrees of divergence in higher dimensions can be higher. For example, 
if a four- dimensional integral produces logarithmic divergence, the corresponding five- 
dimensional divergence is expected to be linear divergence for a naive cutoff regularization. 
This linear divergence can be missed in a dimensional regularization. It needs to be 
examined whether such a regularization dependence occurs for physical effects. Finally 
for brane terms, it was found that loop effects of bulk fields produce infinite contributions 
to require renormalization by couplings on branes [I]. Brane terms can be mass and 
kinetic energy terms and higher derivative operators can be needed as counterterms for 
loop corrections [3]-[8]. Research on effects of brane terms on mode functions and mass 
spectrum has also been developed in the literature [9]- [15]. In principle, the coefficients of 
brane terms and higher- dimensional operators can be energy-scale dependent. In addition 
to divergence, the finite part of quantum corrections needs to be extracted for physical 
quantities. 

In examining quantum loop effects, two-point functions include nontrivial information 
in higher- dimensional field theory. In 4 theory in flat five dimensions on the S l /Z 2 , 
at one-loop level, there is no wave function renormalization. At two-loop level, diver- 
gences for wave function renormalization appear not only for p 2 but also for (p 2 ) 2 . This 
higher-derivative term gives important information for consistency of the theory. The pre- 
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dictability of the theory with (p 2 ) 2 term requires the ultraviolet cutoff orders of magnitude 
larger compared to the compactification scale [16J. One-loop wave function renormaliza- 
tion appears in Yukawa theory. On an orbifold, fermion self-energy has divergence on 
the branes pQ. The corresponding fermion kinetic term not only in the bulk but also on 
the branes need to be included in the starting action integral. At first sight, the relative 
coefficient for the bulk and brane kinetic terms seems arbitrary. It is the case at tree 
level of the action integral. However, for the energy-dependence of physical effects, it is 
crucial to include quantum effects. When one extracts such effects, it must be treated 
carefully whether the relative coefficient is only apparent ambiguity. An analogy lies in 
four-dimensional renormalizable field theory where coupling constants in tree-level action 
integrals are free parameters. After renormalization, they are physical and momentum 
dependent. The difference is that in higher- dimensional field theory full information of 
non-renormalizable interactions is unknown. To proceed phenomenological study without 
such an information, one might be content to assume that the relative coefficient can have 
various values. In this sense, the coefficient is akin to a free parameter. One important 
point is to pay attention to physical consequences. It is nontrivial whether all the non- 
renormalizable interactions are necessary for making a physical prediction at low energies. 
Then one needs the relative coefficient that is not a free parameter and its behavior that is 
found at low energies, whereas the momentum-scale dependence of the relative coefficient 
has not been examined even in the leading level in the literature. 

In this paper, we examine the momentum-scale dependence of fermion self-energy as- 
sociated with wave function renormalization in Yukawa theory in flat five dimensions on 
the S l /Z%. At one-loop level, wave function renormalization is needed as in four dimen- 
sions. One-loop divergence can be subtracted with only two renormalization constants in 
the bulk and on the branes. We show that the bulk and brane parts of the self-energies 
are uniquely determined with requiring physical conditions. Our physical condition is 
that fields with their physical masses obey usual Feynman propagators. This is sufficient 
because a relative difference between left- and right-handed components for fermions is 
related to the difference between the coefficients of the bulk and brane kinetic terms. With 
this procedure, the momentum-scale dependence of the renormalized self-energies is given 
and the distribution of the bulk and brane parts can be found at low and high energies. 
While the loop-momentum integral depends on the regularization scheme, the renormal- 
ized self-energy is scheme-independent. We find that the ratio of the bulk contribution 
to the brane contribution depends on the momentum-scale in a very mild way, so that 
the relative coefficient of bulk and brane kinetic terms can be regarded as approximately 
constant for the leading quantum effect. 

The physical condition given here is applicable to remove ambiguity in various orbifold 
models. Following the recent discovery of marginal and interacting operators in models 
with extra dimensions [17] , we also give the low-energy behavior of bulk and brane kinetic 
terms by the distance-rescaling for integrating out the shell of high-momentum degrees 
of freedom. It would be important to speculate what arises at higher-loop level. The 
regularization-scheme independence may not be kept beyond the leading level. We discuss 
a viewpoint of the regularization-scheme dependence for physical quantities at higher-loop 
level. 

The paper is organized as follows. In Section [21 the model with brane kinetic terms 
is given. The basic idea of our proposal to determine the relative coefficient is described. 
In Section [31 the one-loop divergence is given and is subtracted with the corresponding 
counterterms. In Section HI the physical condition is given for renormalization. The 
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momentum dependence of fermion self-energy is examined. In Section El we give the low- 
energy behavior of bulk and brane operators by the distance-rescaling. We conclude in 
Section [6] with some remarks. A five-dimensional Yukawa theory on the orbifold S 1 /Z 2 is 
given in Appendix [A] The method we employ for calculation of quantum loop corrections 
is exemplified in Appendix [B] At one-loop level, bulk and brane divergences in the five- 
dimensional Yukawa theory are found. Details of mode functions and their orthogonality 
and normalization are shown in Appendix O 



2 Model: brane kinetic terms and mode functions 

We consider a theoretical improvement of a five- dimensional Yukawa theory on the orbifold 
S 1 /Z 2 . The notation is given in Appendix |A] The theory has divergence on the branes. 
The corresponding counterterms are needed. We focus on the effect of brane terms on 
the kinetic energy term. Brane kinetic terms need to be included in the beginning, 

dy [A${i$)1> + B$ L {$)il> L (8(y) + 5(y - L))] , (2.1) 

where the integral | lim e _ i> o f-il e dy is denoted as J dy. The factors A and B are unknown 
constants. Either of these constants, for example, A can be deleted by redefinition of ip. 
Then the equation (12.11) reduces to 

dy + «o M0)i>L (<%) + S(y - L))] . (2.2) 

The coefficient ao should become momentum-dependent after a renormalization. The 
directly-related radiative effect is fermion self-energy. This is given by the sum of bulk 
and brane contributions, 

J VipVipV^ipipe iId4xC , (2.3) 

where the Lagrangian C consists of bulk and brane terms and it is denoted as the four- 
dimensional effective Lagrangian with a Kaluza-Klein decomposition. The basic idea to 
determine the coefficient a is to require that the propagator of the fermion with a physical 
mass M at a renormalization point is given by 

Pl + 7 u , ■ P Ri ( 2 - 4 ) 



j>- M + ie -p-M + ie j> - M + ie 

where Pl and Pr denote projection matrices for left- and right-chiralities, respectively. In 
the equation (J23J), left and right contributions are required to be equal at the renormal- 
ization point. As for the Lagrangian, left and right components have different terms as 
in Eq. (12. 2 j) and the resulting radiative corrections are expected to be different between 
left- and right-chiralities. It is the relative coefficient a that removes this difference by 
the renormalization, while the common part of divergence is removed by the bulk renor- 
malization constant. We will show this occurs in the following sections. 

To treat the renormalization, we define the wave function renormalization factor and 
the rescaled field, 

tj) = Z 1/2 tfj r . (2.5) 
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Substituting the rescaled field into the Lagrangian (12. 2 p yields 

dy [ZTp r (i$)ip r + a Z4i rL (i$)tp rL (5(y) + 5(y - £,))] 
dy [ip r (i$)?p r + a^ rL {i(fi)i\) rL (S(y) + 5(y - L)) 

+S Z A{0)A + WrLWWrL (%) + Kv ~ L))] , (2.6) 

where a is a renormalized brane coupling and the renormalization constants are denoted 
as Sz = Z — 1 and 5 a = a^Z — a. Hereafter we will omit the subscript r to express the 
rescaled field. 

We consider behavior of the coupling a and the self-energy by loop effects for the 
fermion in the Lagrangian 

C = J dy [#7^ + *7 5 W + I (d^ f - l - («9 5 0) 2 - #V0 

+aMi^L(6(y) + S(y-L))], (2.7) 

and the counterterms. The equations of motion for the fermions are 

d 5 iP L + ia»d^ R = 0, (2.8) 
i^d^ L - d^ R + a ia^d^ L ■ (%) + S(y — L)) = 0. (2.9) 

The mode expansion is given by 

oo oo 

Mx,y) = ^2f n (y)^Ln(x), ipii(x,y) = ^ g n (y)ip Rn (x) . (2.10) 

n=0 ra=l 

The Dirac equations icr^d^ipun^x) = m n ipLn{.x) and io^d^Lnix) = m n ip Rn (x) are fulfilled 
by the four- dimensional fields. The two equations have the identical m n for n > 1. Details 
of derivation of the mode functions are given in Appendix The orthogonality is given 
by 

L 

dy fn(y)f m (y) + ^ (/n(o)/ m (o) + f n (L)f m (L)) = 5 nm , (2.11) 

L 

dy 9n(y)g m (y) = S nm . (2.12) 

From these equations, the Lagrangian for the fermion is written in terms of the four- 
dimensional fields as 



L+e 

dy [$(zY*0„ + *7 5 W + aMtf)M5(y) + % - £))] 

L+e 



x[ia -Oxl + YI Whn ia ■ d ^Ln + ipjiJcr • dip Rn - m n ip Rn ijj Ln - m n tp[ n ^ Rn )(2. 13) 



n=l 



The kinetic terms are diagonal with respect to Kaluza-Klein modes due to the orthogo- 
nalities of f n and g n . 
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For the simplest five- dimensional Yukawa theory divergence for non-diagonal compo- 
nents with respect to Kaluza-Klein modes are radiatively generated, starting from diagonal 
kinetic terms. Explicit equations are given in Appendix [B] The equation ( 12. 13[) is diagonal 
for modes. Also in the case with brane kinetic terms, radiative corrections are expected 
to give rise to non-diagonal components. To subtract this divergence, the Lagrangian 
terms in the mode expansion for counterterms need to have non-diagonal components 
with respect to Kaluza-Klein modes. By the mode expansion, the equation (12.61) is 

J^L^ ■ dip Rn + tp[ n ia ■ d*fj Ln ) 
n 

n 

+ E " a5 ^Lm^ ■ dlpLn (/n>(0)/ n (0) + f m (L)f n (L)) , (2.14) 

m n 

where ipw = X- F° r the rescaled fields (the subscript r has been omitted), the kinetic 
terms are diagonal with respect to Kaluza-Klein modes. In the last line, the counterterms 
have off-diagonal components. They are nonzero as we will see below. 

For one-loop calculation, the sum over mass for internal lines must be performed. This 
is difficult when the mass eigenvalue is given in terms of the mass quantization condition 
in the form of a function. Focusing on identifying effects of a, we treat the fermion mass 
and scalar mass at the first order of a as 

mi /_ a\ M , 2 . A nir 



m n = T {l-j)+0(a i ), < = T i (2-15) 



respectively. The corresponding mode functions are given by 



™ -M 1 -^) - + ™ (i - S 8111 } + i2M) 



for < y < L. From these mode functions, the brane counterterm in Eq. (12.141) is 
proportional to 



4 / \ I + ( — I ) n+m 

/»(0)/ m (0) + fn(L)f m (L) = - (l - \ ) + 0(a 2 ). (2.19) 

The off-diagonal components are nonvanishing when the sum (or the difference) of n and 
m is an even number. 

In the present model, the only interaction is the Yukawa interaction. The Yukawa 
coupling for the interaction (i/jLn'4 } Rm4>e + ^Rn^Lm^i) in the four- dimensional Lagrangian 
is 

(-9) I dy fn(y)g m (y)g}(y) = —J== (i - j-) (Se, m+n + $e+n, m - h +m , n ) 



-g) a(n — m) 
L 



I _|_ f ^y+n—m _|_ / ]\l—n+m 



+ 



_ 2(£ + n - m) 2(£ - n + m) 
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£ + n 7^ m £ + m 7^ n 



—q) a(n + m) 



2L 

2\ 



+ 



-r 



_ 2(£ + n + m 



2(£ — n — m) 
n + m ^ £ 



+0{a 



(2.20) 



where £ + n ^ m indicates that the term [l + (— \y +n m J / [2(£ + n — m)\ does not exist 
for £ + n = m. The other indications £ + m 7^ n and n + m ^ £ are similar. 



3 One loop divergence and subtraction 

For the model given in the previous section, we calculate fermion self-energy. The method 
is given in Appendix [B] We discuss correspondences between the divergences and the 
counterterms in the Lagrangian (12.141) . 



3.1 Divergent part 

The first diagram we calculate is the one to have the left-handed fermions with mode n 
in the external lines. The diagrams are shown in Figure [U The internal lines are taken 

s k+n s 

n ; s+n n n ■ k n n n-s '• n 

(a) (b) (c) n-s>0 



s+n n+2s 



s+n 

s n+2s 



(d) 



(e) 



n+2s 



s+n 
s 



(f) 



n+2s 1 s+n 



(9) 



Figure 1: Self-energy diagrams for ipi. 



as ipRm and fa. The modes m and £ are summed. With the equation (12.20j) . the self- 
energy up to the propagators for the external lines is calculated similarly to the one for 
the diagrams in Figures [5] and [TJ The divergence of j) term is found as 



2L(4tt^ 



(1 - 1) at) 



" - log(A 2 i 2 ) 

7T 



(3.1) 



In the limit a — » 0, this reduces to the simplest-model result flB.13j) for ((a) + (b) + (c)) in 
Figure [TJ The liner divergence AL and the logarithmic divergence log(A 2 L 2 ) correspond 
to bulk and brane terms, respectively. This is similar to Eq. (IB. 18j) in the simplest model. 



6 



The divergence (13.11) is subtracted with a linear combination 5z and 5 a for the diagonal 
component in Eq. (I2.14p . For the case with external right-handed fermions, divergent 
part includes a similarly to Eq. (13. ip . For both of left- and right-handed fermions, the 
bulk divergence needs to be subtracted with the identical counterterm. 

The next diagram is the one with ipin arid ipLs in the external lines and with ipRm and 
(pi in the internal lines where the modes m and I are summed. In this case, the divergent 
part for j) term is 



g 2 1 + (-!)"+« 

m n 



2L(4vr 



3a \ , , , 2 2 4a AL 



- 1 - — log(A 2 L 



Pl (3.2) 



2LJ toV ' 3Lv^F_ 

The linear divergence AL appears as effects of non-zero a because the Yukawa coupling 
(I2.20p violates Kaluza-Klein number conservation at tree level. To subtract this diver- 
gence, the bulk kinetic terms in the original action integral seems to need non-diagonal 
components with respect to Kaluza-Klein modes. Introduction of the bulk off-diagonal 
components might be accomplished when the five- dimensional action integral is regarded 
as an effective action in even higher dimensional theory. While such a modification is 
possible, it would be important to identify how the bulk and brane kinetic terms in an 
effective theory vary their leading contributions with the change of momentum scale. At 
least, it is found that the divergences in the limit a — > are completely subtracted with 
the counterterms in Eq. (I2.14p . We leave a further modification of the original action for 
future work. 



3.2 Finite part 

So far we have focused on the divergent part in radiative corrections. To obtain momentum 
dependence of wave function renormalizations, the finite part need to be derived. 

We consider the sum ((a) + (b) + (c)) in Figure HJ From Y^L _oo term, we find 
linear divergence in the bulk with a cutoff regularization. This is the five-dimensional 
correspondent of a four- dimensional logarithmic divergence. On the other hand, the same 
four- dimensional divergence is expressed as r(2 — d/2)|d =4 _ e = T(e) ~ 1/e in a dimensional 
regularization. In five dimensions, T(2 — d/2)\ d=5 = T(— 1/2) = —l^pH is finite. The 
momentum integral seems regularization-dependent. If physical quantities are sensitive 
to the way of the regularization, the predictability would be lost. In the present case, 
it will be found that regularization-dependent divergent terms are independent of the 
momentum scale up to the overall f. Such a momentum-independent shift is unphysical 
for the renormalization. In Eq. (IB.12p . the sum part is written as 




+ lp 2 =0j 



(3.3) 



Only the last term diverge and it corresponds to a constant shift. As seen in Eq. ( IB. 181) 



all the divergence for this part and the others can be subtracted with only the two factors 
5z and 5 a . In the contribution ((a) + (b) + (c)), the momentum-dependent part is obtained 
as 

,2 /•! 



'if) 



2L(4vr 



|2 



[ dxx -log(A 0+ L 2 )-log(A n+ L 2 )+ V logf A A f + ) 



Pl- (3.4) 
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For the external ip n and ip n +2s, the self-energy ((d) + (e)) is 



(if I dx \og{A(m n+s ,m s )L 2 ) ■ P L , 
Jo 



(3.5) 



where only the momentum-dependent part has been given. The self-energy ((f) + (g)) has 
the same value. For the right-handed external lines, the momentum-dependent part is 



fr 



2L(4vr) 2 



„1 oo 

if j dxx log 

Jo „_ _ 



A, 



A 



s+\p 2 =0 



■Pr- 



(3.6) 



In the next section, we will examine momentum dependence of wave function renormal- 
izations for bulk and brane terms at the leading level. 



4 Momentum dependence of self-energy 

Now we perform the renormalization. The full propagator G n ^ is related to the one- 
particle irreducible self-energy £ nj £ and the tree propagator S n n as 



(4.1) 



where £ 



E j>n and G r 



G 



Ge tn . The solution at the leading level is given by 

Sn.n 



n,n 



1 — <? y 

1 u n,n'—"n'i 



G n / — S ntn Ti n ^Se } £ 



(4.2) 



The diagonal and off-diagonal parts are separated. From the Lagrangian (I2.14p . the 
renormalized self-energy — zS, the one-loop self-energy — iEW and the counterterms are 
related to each other as follows: 



i'- 1 !. n.n „ 

-iE Rn>n (fP R 

-iS Ln+ 2s,n(^)-PL 
-iE L „ n+ 2 s (f)P L 



Pl = -iE 



lIS1>)Pl + WzPl + iffl a - a5 z )-P Ll 



-z'S 



(i) , 

R n,n V 
(1) 



)P R + ij>8 z P R , 



-^^MPL + Wa-aSz^PL, 



Ln,n+2s(l>)PL + ~ a5 Z )-P L , 



(4.3) 

(4.4) 
(4.5) 

(4.6) 



where the subscripts L and R label left- and right-handed chiralities for external fermion 
ip with leaving a specification for anti-fermion ip, respectively. Physical conditions for the 
self-energies are given in the following. For the propagator for fermion with a physical 
mass to keep the factor in the numerator i, the renormalization conditions for the diagonal 
part are imposed as 



L n,n 



df> 



0. 



R n,n 



i>=m n 



0. 



(4.7) 



j}=m„ 



where m n is the physical mass for the propagator 



j>-m n - t L , 



Pr 



- m n )(l - t'{m n )) + 0{{i> - m n f) P L , (4.8 
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with the same form of equation for S^„ jn . From the condition for in Eq. (14. 7p . the 

counterterm 5z is obtained as 



(i) 

Rn,n 



(4.9) 



=m n =p5 



where p$ is the Fourier transformation of extra-dimensional derivative. The equation for 
the counterterm should be interpreted in terms of p§ rather than the form dependent on 
Kaluza-Klein modes. From this equation and the condition for the left-handed component 
in Eq. (14.71) . the other counterterm is obtained as 



(5 a - a5 2 



,(1) / 
J Rn,n\ 



(4.10) 



j>=m n =p5 



The renormalization conditions for the off-diagonal external lines n + 2s, n and n, n + 2s 
need to be imposed as £r,n+2s,n(P5) = ^Ln,n+2s(p^) — 0, so that the field for the propagator 
(14. 8p is in the mass eigenstate. In these equations, the Kaluza-Klein mass-dependent part 
in ^Ln+2sn an d ^Ln+2sn should be read with the Fourier-transformed derivative p 5 . 

These self-energies can be decomposed into bulk and brane parts. The bulk part is 
defined from the contribution for the right-handed fermion as 



2L(4vr) 



bulk( 



^(1) / 
J Rn,n\ 



(4.11) 



Here the correspondent with hat is defined similarly. From the equations (14. 4p . (14. 9p and 
(14. lip , the renormalized Sbuik is 



J bulkl 



o 



l oo 

doc vC ^^^^^ 



log 1 + 



x(l - x)(pl - p 2 )L 2 
(sit + (1 - x)p 5 L) 2 



+ 



2x(l — x)p'p 5 L 2 
sir + (1 - x)p 5 Lf 



.(4.12) 



Next the brane part is defined from the difference between the contributions for the left- 
handed and right-handed fermions as 



9" 



2L(4vr) 



(1) (4) 

brane \"1 



_ E (D 



Rn,n 



(4.13) 



Here the correspondent with hat is defined similarly. From the equations (14.31) . (14. 4p . 
(I4.10p and (I4.13p . the renormalized Sbrane is 



brane \ 







dx ( log 



-x(l — x)p 2 + (1 — x)p\ 
(1 - x) 2 pl 



+ log 

+ 



— x(l — x)p 2 + (4 — 3x)pl 
(2 - x) 2 pl 



dx 2x 2 (l — x)^p5 







+ 



(l-x) 2 p 2 (2-x) 2 p 2 



(4.14) 



Thus the self-energies (I4.12p and (14.141) have been determined uniquely. 

The bulk and brane parts of the self-energy is numerically estimated. For the equation 
(14. 12p and (I4.14p . the self-energies are singular for x = 1. This is the same as in infrared 
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Figure 2: The momentum dependence of EbuiktjO' where p$ = 1/L and — 10 3 < s < 10 3 . 

divergence for massless photon in the four-dimensional quantum electrodynamics. If the 
five- dimensional scalar field is massive, this infrared divergence is avoided. For the numer- 
ical analysis, we take the values of p$ and the five-dimensional scalar mass as the typical 
dimensional quantity 1/L. By defining Q = (^)L)/(p 5 L), we examine the momentum de- 
pendence of the bulk and brane self-energies. The momentum dependence of Sbuik(^) is 
shown in Figured For modes of the Poisson summation, we take the sum — 10 3 < s < 10 3 
as the contribution for \s\ > 10 3 is much less than 0(%). For Q ~ 0.4, the linear and 
logarithmic parts of Eq. (I4.12p are comparable. For 0.7 < Q < 1, the linear part is dom- 
inant. For Q — 1, i.e., j> = ps, the logarithmic part is vanishing. The renormalized bulk 
self-energy is written as the linear term and the others such as Ebuik ~ 0.7Q + £ b 1 ™j 1 1 [ mear , 
where Sb u ik|Q=i — 0.7. For 10~ 5 < Q < 1, the bulk self-energy Sbuik changes the value by 
about 50%. The momentum dependence of Sbrane(^) is shown in Figure El The behavior 



0.30 



0.25 



0.20 - 




Figure 3: The momentum dependence of Sbrane(^), where the value for p 5 is taken as the 
same as in Figure |2J 

is similar to the case of Sbuik- The brane self-energy is written as the linear term and 
the others such as Sbrane ~ 0.3Q + S£°^ 1 ( f ear , where Sbrane|Q=i — 0.3. In addition to each 
momentum- dependence of the bulk and brane parts, the momentum-dependence of the 
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ratio is important. The ratio of the brane contribution to the bulk contribution is shown in 
Figure HI For 0.7 < Q < 1, the ratio is dominated by the linear term and is approximated 



0.4180 



J 0.4175- 



0.2 0.4 0.6 0.8 1.0 

Q 



Figure 4: The ratio (Sbrane/Sb u ik)(^), where the value for p 5 is taken as the same as in 
Figure [H 

as Sbrane/Sbuik ~ OAQ + (nonlinear), where 0.3/0.7 ~ 0.4. The nonlinear terms have 
logarithmic dependence of Q and the change of the value against the momentum- scale is 
small. For 10~ 5 < Q < 1, each of Sbuik and Sbrane changes the value by about 50%. On 
the other hand, the change of the value for the ratio is less than 0.1%. The momentum 
insensitivity of the ratio is due to the analogous form seen in Sbuik and Ebrane- 

From the result for the self-energy, it follows that the ratio of the brane contribution 
to the bulk contribution is almost constant with respect to the momentum-scale. When 
(5 a — a5z)/Sz is constant, the S a is represented in terms of 5z as 5 a = (a + constant)^. 
At the leading correction, a = 0. Then the Lagrangian (I2.6P can be represented as 

J dy [Vv(^)W + S z {A(0)A + 0A2i> rL (i^ rL (5(y) + 5(y - L))}] , (4.15) 

where 0.42 is the value of the ratio for the region in Figure HI Therefore the relative coef- 
ficient between the bulk and brane kinetic terms has been determined without ambiguity. 

5 Bulk and brane operators in the distance-rescaling 



In this section, we discuss the issue of relevant, irrelevant and marginal operators for bulk 
and brane terms to identify a general low-energy behavior of quantum corrections in the 
field-theoretical context. 

In Ref. [T7], marginal and interacting operators in quantum field theory with extra 
dimensions were discovered for the Randall-Sundrum spacetime whose metric is given 
by [H1CL9] 

ds 2 = \ (ri^dx^dx v - -j^dzA , (5.1) 
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where k is the curvature of the five-dimensional anti-de Sitter space. For the four- 
dimensional rescaling x' = bx where b < 1, the rescaling of z is given by 2' = z/b. 
The action integral 



(Fx 



dz 
kz 



1 1 1 

—ipij^d^ + t^(<V>) 2 - -jgipip<j> 



is rescaled into 



dz' 
kz' 



(5.2) 



(5.3) 



Here the field redefinition is given by if)' = b~ 3 ip and 0' = b^ 2 <p. The coupling constant is 
given by g' = b°g and the corresponding interaction is a marginal operator. 

If the spacetime is flat with the rescaling x' = bx and y' = by, the action integral 



d A xdy 



is rescaled into 



d A xdy 



(5.4) 



(5.5) 



Here the field redefinition is given by ip' = b~ 2 ip and <\>' = b~ 3 ^ 2 (p. The coupling constant 
is given by g' = b x l 2 g and the corresponding operator is an irrelevant operator. 

Now we consider the sum of bulk and brane terms. In the flat spacetime, the action 
integral is given by 



d 4 xdy C hu i k (ip, 4>) + / d 4 xdyC hTane (ij, <p)5(y - y { 



(5.6) 



where yi denote the positions of branes. For the rescaling in the flat spacetime, the action 
integral is written as 



d 4 x'dy'b- 5 £ bu i k (&y , & 3/ V ) + / d 4 x'dy'b~ 5 £ brane (6V, & 3/ V )<%' - v[)h. (5.7) 



The brane terms are multiplied by the factor b compared to the bulk terms. This shows 
that in the flat spacetime brane kinetic terms are small compared to bulk kinetic terms. 

In the Randall-Sundrum spacetime, the action integral corresponding to Eq. (15.61) is 
given by 



d A xdz^J det0Mjv£buik(^, & 9mn) + / d A xdz^J detg MN C hra , ne (tp, <fi; g MN )kz5(z - ^)(5.8) 



For the rescaling in the Randall-Sundrum spacetime, the action integral is written as 



d 4 x'dz'^/detg' MN b 8 C hu i k (b 3 ip' , b 2 <p'; b 2 g' MN ) 
+ / d 4 x'dz'^/detg' MN b- 8 C hr&ne (b 3 i>', b 2 <j) f ; b- 2 g' MN )kz'6(z' - z[). (5.9) 
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Because kzS(z — Z{) = kz'5(z' — z'^b , the brane kinetic terms are the same order as the 
bulk kinetic terms. 

Our method to determine the momentum-dependence of the bulk and brane contri- 
butions has been given in the flat spacetime. In the flat spacetime, the bulk interactions 
and the brane kinetic terms tend to be irrelevant operators. Thus our explicit analysis 
may be regarded as the proposal of an idea in a simplified model. On the other hand, 
in the Randall-Sundrum spacetime, these operators tend to be marginal operators. In 
the four-dimensional theory, marginal operators are renormalizable terms and low-energy 
quantum corrections are consistently treated without an additional ultraviolet comple- 
tion. The idea needs to be examined further in models where the bulk interactions and 
the brane kinetic terms are marginal operators. 

6 Conclusion 

We have examined energy-scale dependence of fermion self-energy associated with wave 
function renormalization in Yukawa theory in flat five dimensions on the S 1 /Z 2 . At 
one-loop level, wave function renormalization is needed as in four dimensions. One-loop 
divergence can be subtracted with only the two renormalization constants 5z and S a in the 
bulk and on the branes. We have shown that the bulk and brane parts of the self-energies 
are uniquely determined with requiring physical conditions. The physical condition is that 
fields with their physical mass obey usual Feynman propagators. In the present setup, for 
the right-handed component the bulk part of the self-energy is immediately fixed. The 
left-handed component determines the other renormalization factor. With this procedure, 
momentum- scale dependence of the renormalized self-energies has been given. While the 
loop integral has linear divergence in the cutoff regularization, the renormalized self-energy 
does not give rise to regularization-scheme dependence. As an explicit equation, the bulk 
part of the renormalized self-energy is rewritten as 

From this equation, it is seen that a constant shift of S^ ulk is unphysical. The case of 
Sbrane is similar. We have also found that the ratio of the bulk contribution to the brane 
contribution depends on the momentum scale within less than 0.1%, so that the relative 
coefficient of bulk and brane kinetic terms can be regarded as approximately constant for 
the leading quantum effect. The value of the relative coefficient has been estimated as 
0.42 for the parameter set given in Section HJ Because the value is not extremely small, 
it would affect the prediction in phenomenological applications such as collider physics. 
The quantitative estimation depends on the model. 

Analysis beyond the leading level is an open question. In the equation f)3.2p . divergence 
for the off-diagonal component with respect to Kaluza-Klein modes has been found. To 
subtract this divergence, the original action integral would need the corresponding coun- 
terterms. As mentioned below Eq. (13.21) . linear divergence for the off-diagonal component 
may correspond to five- dimensional terms in higher- dimensional spacetime. The necessity 
of new local operators order by order is nothing but non-renormalizability. A possible 
way to proceed is to show that such local operators are irrelevant. In the present analysis, 
the relative coefficient of the brane terms to the bulk term is 0.42 which is not negligible. 
The effects of new operators should be treated carefully. As for regularization-scheme 



(6.1) 
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dependence, it needs to be checked whether the scheme independence is kept at higher- 
loop level. Because the loop-integral is scheme dependent even at the leading level, the 
momentum-dependent part for any physical quantities might be contributed from such a 
effect. A viewpoint of this issue is given in the following. When our interest is low energy 
behavior of physical quantities, details of the scheme- dependent part may not be required. 
Then we do not need to show that the prediction is exactly scheme independent. Rather, 
it is sufficient that the scheme-dependence is negligible at low energies. In addition, our 
analysis has focused on wave function renormalization for two-point functions. To define 
a physical mass, mass renormalization should be treated appropriately. 

Finally we emphasize that the physical conditions given here are applicable to remove 
ambiguity in various orbifold models. We have also shown that both of the bulk and 
brane kinetic terms are marginal operators in the Randall-Sundrum spacetime. It would 
be worthwhile to identify the momentum-dependence of bulk and brane contributions in 
such models. 
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A Yukawa theory on an orbifold 

We consider a five-dimensional Yukawa theory on the orbifold S 1 /Z 2 . The bulk spacetime 
is flat with the metric rj MN = diag(l, —1,-1,— 1,-1). The action integral for the Dirac 
fermion ijj(x,y) and the real scalar field <p(x,y) is 



S 



d 4 x ■ - 



dy 



(A.l) 



Greek indices /i run over 0,1,2,3 and fifth index is denoted as y. The gamma matrices are 



a 1 



-Z7 



(A.2) 



where = (I2, a 1 ) and <x M = (— 12, c l ) with the Pauli matrices. Dirac conjugate is given 
by ip = t/^7°. The fields have period 2L with respect to the extra dimension. The orbifold 
boundary conditions for the fields are chosen as 



~y) = h 5 i>(x, y), ip(x, L-y) = i^ipix, L + y), 
<f>(x, -y) = -<P(?, y), <f>(x, L-y) = -<p(x, L + y), 



(A.3) 
(A.4) 



for the extra-dimensional fundamental region < y < L. With this boundary conditions, 
the fermion has mode expansion as 



1 °° [2 

-7= Xl{x) + y 1 ^Ln(x) cos (m n y) 

* n=l 



i>L{x,y) 



ipn(x,y) = Y \ T^Rn( x ) sin ( m ny) 



n=l 



(A.5) 
(A.6) 



with m n = rnr/L. The projection matrices are given by Pl = [I2 — ( — *7 5 )] /2 = 
diag(l, 1,0,0) and Pr = [I2 + (— «7 5 )] /2 = diag(0,0, 1, 1). The mode expansion of the 
scalar <p(x,y) is similar to Eq. flA.6j) . The mode functions and the above trigonomet- 
ric functions are periodic for y — >■ y + 2L. From the properties of orthogonality and 
completeness, quadratic terms in the action integral (IA.1I) are 



dy 



XL 



n=l 



11 

ipnih^d^ - m n )ip n + -(d^n) 2 + -(m n 0„) 2 



(A.7) 



where ip n is composed of ipL n and ipR n as in the left- and right-handed projections for 
ip. The equation ( 1A.7|) is diagonal with respect to Kaluza-Klein modes. The Yukawa 
interaction is written in terms of four-dimensional modes as 



1 f L - °° (-9) 

O / d V (-9)^^ = V —=- [xL^Rn<Pn + 

IJ-L n=1 VL 



1pRnXL<t>n] 



00 00 00 



+ - 1= 1pLn4>Re4>k {Sk+n,£ + <^+n,fc — <W^,n} 



n=l £=1 k=l 

00 00 00 



2L 
-9) j. 



+ -/= IpRnlpLi&k {$n+£,k + 5, 



71 = 1 £=1 fc=l 



k+£,n ~ On+fc 



a . (a.8) 
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The Yukawa interaction involves the exchange between Kaluza-Klein modes. Originally 
the kinetic terms are set only in the bulk. In the next section, it will be shown that 
divergence corresponding to brane kinetic terms is radiatively generated. 

B Method of loop calculations 

As in four dimensions, wave function renormalization is expected to appear for the fermion 
at the one-loop level. The appearance of divergence on the branes in a similar setup is 
shown in Ref. PQ. For the calculation of quantum loop corrections, the steps we perform 
contain the Kaluza-Klein mode expansion, the summation of diagrams, the replacement 
of fractions with a formula for the Gamma function and the representation with the 
Poisson summation for the summation of the Kaluza-Klein modes. In this appendix, we 
give explicit equations at each step in the method to calculate self energy for the action 
integral (lA.ljl . The method will be used further in Section [3j 

The first diagram is the self-energy for the left-handed fermion zero mode. The external 
lines are the left-handed zero mode x an d the internal lines are <fr n and ipR n - The internal 
lines are summed with respect to all n. The diagram is shown in Fig. [5j The self-energy 

p-k 

n < 

P /' k \ P 



n 



Figure 5: Self-energy diagram for the left-handed zero mode, 
up to the propagators for the external lines is 

9 2 ^ f d A k I^ + m n 1 = 9 2 ^ [\ [ dH xj> + m n 

L^J (2n)*k?-m* L {p - kf - ml L ^J J (2tt)* [P - A.] 2 ^' 

with a Feynman parameter x. Here i = k — xp and A n = — x(l — x)p 2 + m 2 n . We will 
focus on p 1 terms, i.e., the terms for the wave function renormalization, neglecting m n in 
the numerator. Then the integrand is symmetric under n -H- — n. The j> term in Eq. (IB. 1[) 
is rewritten as 

9 - V fdx [ *M 2 P L -^r d x[ *M 2 P L , (B.2) 

with the Wick rotation £° = i£° E and t = t E . With use of symmetry under n -H- — n, 
the sum in the first term has been taken for — oo < n < oo which is more convenient 
than treating the sum for 1 < n < oo. For this infinite sum, a formula for the Poisson 
summation will be employed. The second term in Eq. (IB.2|) is independent of internal 
mode n. 

Let us calculate the sum part Y^=-oo- The first step is to rewrite the fraction using 
the integral representation of the Gamma function as 

1 -i roo /*oo o 2 

1 1 / . . . ._ ri>2 j A U / . re2 „n „\„2 i n 7r 



[P E + A n ] 2 r(2) 



POO pOO 2 2 

/ dtte~^ +A ^= dtte-^-^ 1 -^^. (B.3) 
Jo Jo 
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The second step is to change the mode index with the Poisson summation, 



e 



L 



(B.4) 



Then the sum part, the first term in Eq. ( IB. 21) . is written as 



2 5 7T§ 



pi x poo 

/ dxx(ij)) / dtr 3 *e x{l - 



-x)p 2 t- 



(B.5) 



where the Gaussian integral for £e have been employed. The next step is to evaluate 
each part for n p = and n p ^ 0. The n p = part includes divergence for t — > 0. The 
divergence of the i integral is evaluated as 



dtt~*(l + x(l -x)p 2 t + ---) 



f =00 
t=A- 



2A + finite, 



(B.6) 



at the cutoff regularization. This momentum integral depends on the regular izat ion 
scheme. The scheme-independence of the renormalized quantity will be shown in Sec- 
tion HI The divergence of the n p = part is 



2 5 7T 



dx x{if) ■ 2 A P L 



2 JO 



9 2 
2L(4vr) 2 



(B.7) 



From the mass dimensions for the fields and coupling [ip] = [mass] 2 , [0] = [mass] 3 / 2 , [g] = 
[mass] _1//2 , the bulk divergence for the wave function is expected to have the coefficient 
g 2 A where [g 2 A] = [mass] and it is consistent with Eq. (1B.7|) . The n p ^ part is 



2 3 7T 



2 JO 



1 °° I— 

dxx{if) ^— e - 2n ^- x ^ x)p2 P L 
2n p L 



(B.8) 



where the integral for t has been performed with help of expressions for Bessel function 
K v . If a function F{x) is symmetric under the exchange x -H- 1 — x, F(x) = F(l — x), 
then Jp 1 dxxF(x) = L dx (1 — x)F(x) = (1/2) dx F(x). The equation (1B.8|) becomes 



2L(4tt)^ 



dx {ifi) log 



1 -e 



-2Ly/ -x{l-x)p 2 



■Pi 



(B.9) 



where the sum for n p has been performed with the Taylor series expansion log(l — X) 
Now we move on to the second line in Eq. ( IB. 21) . 



(A part) = - 



2L 



dx 



d 4 £E x(i$) 

(2tt) 4 \e E + a ] 



2 PL, 



(B.10) 



where A = —x(l — x)p 2 . This is a usual four-dimensional integral. For the cutoff 
regularization, Eq. (IB.lOj) is given by 



2L(4vr 



i2 [ij>) ( log(AL) - i \og{-p 2 L 2 ) + 1)'Pl- 
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(B.ll) 



This divergence is expected from the dimensional counting [g 2 ] = [mass] -1 and [1/L] = 
[mass] . For the self-energy with the field x i n the external lines, we have found linear diver- 
gence in the summation part Yl^-oo anc ^ logarithmic divergence in the non-summation 
part. 

The self-energy with the field i/jl in the external lines is analogous. Here there are 
various Kaluza-Klein modes to contribute. The internal lines are composed of and 
ijjR. The Kaluza-Klein modes for the external and internal lines are labeled as shown in 
Figure [TJ The momentum flow is the same as in Figure |5j With the same procedure as 
in obtaining Eq. (1B.2j) from Eq. (IB.lj) for the diagram in Figure |5j the part j> in the sum 
of the diagrams (a), (b) and (c) in Figured] is obtained as 



n 2 00 r 1 r 

((a) + (b) + (c))|, = £ J dxj 



d A 



00 

2 /•! 



~k 1 (lr 



'0 

-.2 rl 



(27T) 4 [£2 - As+ 

d A i xj) 



2L 



dx 



(27T) 4 [P - A(H 

d 4 £ xj) 



(27T) 4 [P - A r 



;b.i2) 



Here A s+ = —x(l — x)p + (s + (1 — x)n) 2 7r 2 /L 2 + x(l — x)n 2 n 2 /L 2 . The other A 0+ and 
A )t+ are obtained from A s + as A + = A s+ | s=0 and A n+ = A s+ | s=n . The divergent part 
is given by 



((a) + (b) + (c))| div 



2L(4tt) 



^ - log(A 2 L 2 ) 

71 



■Pi, 



(B.13) 



For the diagrams (d), (e), (f) and (g) in Figure[TJ the external lines have n and n + 2s. The 
modes n and s in the internal lines are fixed for given external lines. These correspond 
to the non-sum part in Eq. (IB. 21) . These diagrams satisfy (g) = (d) and (f) = (e). The 
divergent part is 



(d)| 



div 



div 



4L(4tt) 2 
(f)|div = (e)|div 



if) log(A 2 L 2 ) • P L 



(B.14) 



The self-energy with the fields \ and ipL in the external lines is similarly calculated. 

The other diagrams are the ones that have the right-handed field ippt in the external 
lines. Similarly to diagrams with ip^ in an external line, the Kaluza-Klein modes are 
labeled. For (a), (b), ■ • •, (g), there are the corresponding diagrams (a'), (b'), ■ ■ •, (g'). 
The amplitudes (a'), (b'), (c'), (d') and (g') are given by (a), (b), (c), (d) and (g) with Pl 
replaced by Pr. The amplitudes (e') and (f ) are given by (e) and (f) with Pl replaced 
by Pr and with the change of the overall sign. In addition to these diagrams, there is 
a diagram with ipR in the external lines and with x an d 4> in the internal lines. This 
contributes to the divergent part as 



2L(4tt) 



(#) log(A^) ■ P R 



(B.15) 



The equation (1B.15P and the last term in Eq. (IB. 13[) with the replacement Pi by Pr cancel 
each other. The contributions (d') and (e') cancel each other and (f) and (g') cancel each 
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other. In summary for the right-handed external lines, the remaining divergence is 



{if)—.p R . 



(B.16) 



The logarithmic divergence for the external ipR is canceled out. 

The final step in the loop calculation is to identify Lagrangian terms corresponding to 
the divergence. For all the divergences (1B.13j) . ( 1B.14|) and (1B.16j) . the Lagrangian terms 
are given by 



-i)g 2 al 



2L(4vr) 



i>)XL +^2^n{if)^n 



\0g{kL){x L {if)XL + 2j2^n, L (¥)^n,L 



2L(4tt) 2 

+ 2 ^2 ^2 {$n,L{ifi)i>n+2s,L + ^n+2s,L{if)^n,L) 
n s 

+ ^5^ {Xh(if)^2n,L + 4>2n,L(ifi)XL) 



(B.17) 



From this equation, the Lagrangian terms with four-derivative are written in terms of the 
five-dimensional field as 



— % — — ■ - lim 

2L(4vr) 2 2 e^o 



L+e 



dy 



L+e 



AL 



- log(AL) Mi^ L {5{y) + 5{y - L)) 



(B.18) 



where lim e ^ j_T +t dy G(y) [5(y) + 8{y - L)\ = G(0) + G(L). The diver gence in the brane 
terms for one-loop wave function occurs in the equal size at y — and at y = L, in 
agreement with Ref. [TJ. For the boundary condition ( 1A.3I) . the right-handed component 
ipR does not have the brane term. In the equation ( IB. 181) . the brane divergence includes 
non-diagonal components with respect to Kaluza-Klein modes and the bulk divergence is 
diagonal. 



C Mode functions 

According to the approach given in Ref. 0, we derive f n and g n . The fermion mode 
functions appearing in the context, f n and g n , are f n m and g n m given in this appendix. 
From the equations of motion (12.81) and (12. 9p and Eq. (I2.10p . the mode functions obey 

dyfn = -m n g n , (C.l) 
d y g n = m n (l + a (5(y) + 5(y - L))) f n , (C.2) 

for n 1. From these equations, the second-order differential equation for /„ is obtained 

as 

(d 2 y + ml + am 2 n (8(y) + 8(y - L))) f n = 0. (C.3) 

Unless confusion arises, we will omit the subscript for the Kaluza-Klein mode number n 
for f n . The equation ( 1C.3I) includes two delta functions. It is convenient to write the 
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Figure 6: Regions. Allows depict 5 functions 



function by defining separate regions in a period 2L. The regions can be classified as 
shown in Figure [61 In each region, the function / where is written as 

fi{y) = Ae imv + Be~ imy , (C.4) 

f a ( y ) = Ce imy + De~ imy , (C.5) 

f m (y) = Ee imy + Fe~ imy , (C.6) 

f w {y) = Ge imy + He~ imy , (C.7) 

where A, B, C, D, E, F, G, H are determined by 8 conditions. They are periodicity between 

the regions I and III and between the regions II and IV, continuity at y = —L, 0, L, 
matching of the first derivative at y = —L, 0, L and the normalization. Among the above 
9 conditions, one condition is automatically satisfied if the other conditions are satisfied. 
Now we fix these constants. The periodicity requires 

Ae im{y-2L) + Be -m{y-2L) = gjmy + p^i-my ^ ( C _g) 

C jm{y-2L) + De -m{y-2L) = Q( J,my + He -imy g) 

The continuity requires 

Ae~ imL + Be imL = Ce~ imL + De imL , (CIO) 

C + D = E + F, (C.ll) 

Ee imL + p e -imL = Q^mL + He ~imL ? ^ 

at y = —L, 0, L, respectively. The matching of first derivative requires 

[dyfnfi + am 2 n fnl\y=-L = 0, (C.13) 

[dyfn)n + am 2 n f n ii\ y= o = 0, (C.14) 

[dyfnfui + am 2 n f nm \ y=L = 0, (C.15) 



at y = —L, 0, L, respectively. Here f n i in the second term in Eq. (1C.13j) can be replaced 



by fnii because of the continuity (1C.10I) . The equations (1C.14j) and (1C.15I) are understood 
similarly. In summary, the conditions up to normalization are 

Ae im(y-2L) + Be -im(y-2L) = pjmy + p^imy ^ ( for 0+ < y < (C.16) 

Ce im(y-2L) + De ~im(y-2L) = Qjmy + He ~im V ^ ( for L+ < y < 2 L_), (C.17) 

Ae -imL + Be imL = Ce ~imL + D jmL^ (Q^g) 

C + D = E + F, (C.19) 
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Ee imL + Fe -imL = QjmL + He ~imL^ 

Ce -imL _ De tmL _ ^ + iam y^mL + £^ _ iam yimL = ^ 21 ) 

E — F - C(l + zam) + D(l - iam) = 0, (C.22) 
Ge imi - He~ imL - E(l + iam)e imL + F(l - iam)e- imL = 0, (C.23) 

where 0+ = lim e ^ + e and L T = lim^o-^ T €■ Using the first three equations (1C.16I) . 
(IC.17P and (|C.18p . the fifth condition (IC20[) is automatically satisfied. Hence, the above 7 
conditions and the normalization fix the 8 constants A, B, ■ • ■ , H. Up to the normalization, 
the solution of the mode function / is given by 

fry) = A(e imy + w 3 e- imy ), (C.24) 

f Jl (y) = A(we im v + w 2 e- imy ), (C.25) 

f m (y) = A(w 2 e im y + we^), (C.26) 

f w (y) = A(w 3 e im y + e- my ), (C.27) 

with the mass quantization condition 

civn 

-tanmL = (C.28) 

l-(2f) 2 

which is also written as e~ 2linL = {2+iam) 2 / (2 — iam) 2 . Here w n = (2+iam n )/(2—iam n ). 

The orthogonality of the mode function can be found as follows. The mode function 
satisfies the second-order differential equation ( 1C.3j) . From this equation, the following 
equation is derived: 

= \ lim dy{f m (d 2 y + m 2 n + am 2 n (6(y) + 5{y - L))) f n 

~fn {d 2 y +m 2 m + am 2 m (<%) + <% - L))) f m ). (C.29) 
It is found that the mode function satisfies 

lim(/ m ^/ n ) IVl i +e = A m A n (w m + l)~ immL {w n - l)im n e- im " L 

e— >0 

= ]im(f m d y f n )n,-L+e, (C.30) 

and 

r-0 r rL 



dy fun(y)fum(y) = d(-y) fii n (-y)fn m (-y) = dyf mn (y)f mm (y). (C.31) 

J L JO 



Using these equations, Eq. (1C.29I) becomes 

= (m 2 n - m 2 J y dyf n (y)f m (y) + ~ (/»(0)/ m (0) + /„(L)/ m (L))} , (C.32) 

where lim e _j. Ji +e dy (smooth function) = 0. Therefore the orthogonality is given in 
Eq. ( 12. lip . Here A is normalized as 

2A 2 w 3 (l+ - 4Q 2 ) =1. (C.33) 
V 4 + a?m 2 
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The completeness corresponds to 



fm(y)fm(y') = 2% -y')-a (%) + 5(y - L)) fm(y)f m (y')- (CM) 

m m 

The mode function g(y) is derived from the relation g = —(l/m)d y f. They are 

gi {y) = -At(e m y-w 3 e- m y), (C.35) 

gu (y) = -Ai(we imy - w 2 e~ imy ), (C.36) 

gm{y) = -Ai{w 2 e tmy - we~ my ), (C.37) 

m (y) = -Ai(w 3 e imy - e~ imy ). (C.38) 

The function g has the orthogonality given in Eq. (12.121) . Here + in ± has been adopted 
as w n = e~ imy in ±e~ tmy . In the limit a = 0, these mode functions reduce to the mode 
functions in Eqs. (1A.5I) and (1A.6I) . 
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